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иыхъ дифферентальиыхъ уравнен1й, 


4. И. Брылова. 


ПРЕДИСЛОВЕЕ. 


Эта статья служила предметомъ сообщешя въ союзЪ Морекихь Ин- 
женеровъ на годовомъ его собран! 1917 года и войдеть въ „Извьейя 
союза Морскихь Инженеровь“ за этоть годь, 

Но „Извъеня 0. М. И.“ имБють раепространене лить въ сравни- 
тельно тЪеномъ кругу инженеровъ этой спещальности, вопроеъ-изе объ 
нитегрировани  обыкновенныхь дифференшельныхь уравнен есть 
вопросъ обный, предотавяяющ значительный интересь не только для 
инженеровъь другихь спецзальностей, нп и для физиковь. 

Въ самомъ дьлЬ, къ обыкновеннымь дифференщеальнымь уравие- 
нямъ приводять, всЪ вопросы динамики, когда надо опредълить дви- 
жеше матеральной точки, системы точекъ, твердаго тЬла или епстемы 
такихь тьль. Къ таклмъ же уравненямъ приводять вопросы о равно- 
ве! гибкихь нитей. о равновфеш тонкихь упругихь стержней, о 
формЪ калель жидкости и вообще объ опредъленш поверхности жнид- 
кости подь вщящемь капиллярныхь силь, когда эта поверхвость есть 
поверхность вращеня, вопросы о взаимной индукц и распространени 
токовъ по линейнымь проводнихамь м. т. д. КромЪ этого интегриро- 
ван!е многихь уравнеюй въ частныхь производныхь встрычающихся 
въ математической физикЪ приводится къ пнтегрированно уравненй! 
обыкновенныхф. 

Въ весьма многихь случаяхъ эти уравнешя не подходять поль 
ТБ немноме „типы“, которые интегрируются въ конечномъь видЪ пли 
квадратурахь, и часто бываеть. что для дальнЪйнтаго рЪщевя эгра- 
ничиваются частными случаями вопроса, нногда-же въ самомь урав- 
ненит отбрасывають нфкоторые члены или-же прикимають нъкоторые 
параметры. входяцие въ уравнене, за постоянные, п такимъ образомъ 
завъдомо вводять погрышностн, оцфнить вйяве которыхъ невозможно. 

Воть въ ятихь то случаяхь приближенное численное питегриро- 
ваще уравнен!: притомъ такое. въ готоромь не приходится дифферен- 
нировать входявыя въ спставъ уравненя извфетныя функции перем%н- 
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ной независимой; и предетавляеть важное подспорье при физичеекихь 
изслБдованяхь. Эдфсь эти заланныя функши могуть быть уволимы 
Въ уравнене, пользуясь графическою кривою ихь представлякицею или 
зупирнчески составленною таблицею ихъ значен!. Такимъ образомъ 
уравнене иожеть бытъ взято зъ вид соотв уствующемь опыту; чие- 
деяное рышене его можеть быть выполнено съ точностью, которая 
также соотвЪтетвуеть точности опытовъ, п елЪдовательно, теоретичесше 
и наблюденные результаты явятся вполнЪ между собою еравнимыми 
и позволять съ увфренностью судить о правильности или опгабочноети 
тьхь гипотезъ, которыя послужили для составлешя самого дифферен- 
цальнаго уравненя. 

Методы танего чиеленнаго интегрирован я почти не излагаются 
въ вовременныхь куреахь анализа. пресльдующикь не столько прак- 
личестя прикладныя цёли, сколь отвлеченно математическую. такъ 
сказаль, философекую сторону науки, потому эти методы не пользуются 
общею извЪстностью, и зачастую у учащихся остается оть куреа диф- 
ференщальныхь уравненй! то впечатлЪьве, что если предложенное 
уравнен:е ни въ одному изъ этиповъ“ пе подходить, то еъ нимъ ни- 
чего и подълать нельзя. 

Въ виду этлхь соображенй въ отвЪть на леетное для меня пред- 
дожене проф. Академика П. Ц. Лазарева, принять учаете въ „Архив 
Физическихь Наукъ“ я предложаль помъетить настоящую статью, 
которая по своему содержаннь хотя н не относится къ физик, но 
можеть облегчить приложешя математикн къ изельдованио физиче- 
скихъ вопроеовъ. 

Самую статью я оставить въ томъ видь, въ какомъ она была 
первоначально написана и лишь въ концЪ добавиль одинъ примьрь 
чисто физнческаго уравнензя, относящагося къ опредъленю формы 
капли. 

о поводу этого примфра необходимо сказать, что въ сочинении: 
Ап ацешре {0 1е3Ё 4пе Леопез оЁ саргу асйов.. 5у Ву. Ваз/ов 
УИ ап ехапаЦон оЁ 1Ве шефо4 оЁ ицестамоп ешрюуе@, БУГ. 6. 
„АЧать, Саш 4ке 1883, этоть знаменитый астрономъ излагаеть методу 
интегрировашя лифференщальныхь уравнени! пользуясь конечными 
разкостями. Но надо ныфть въ виду, что Адамоъ быль етоль же зам 
чательный вычислитель какъ п теоретикъ, Объ этомъ свильтельствують 
напримфрь, его стать: „Тье сйешацов оГ Вегпо ав полвегв опа 
В» № Вы пли „№юз оп Зе узше оГ Ешегз сопзгат“. Достаточяо 
хказять, что чиелитель наименьшаго изъ чисель Бернулли имъ вычя- 
сленныхь В„ есть нЪлое число въ 12 цифры. а послёдеяго В» въ 
319 цифрь, что Эвлерона постоянная имъ вычислена сЪ 273 знаками, 
чтобы вилъть, что человькь не обладающёй спезельною способностью 
и необыкновенныхь навыкомъ такихь’вычислевй выполнить не можеть 

„ 
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Возможно, что въ силу привычки къ вычисленямъ еъ весьма 
больою точностью Адамеь и въ излагаемомъ имъ епособъ прибли- 
женнаго интегрировашя стреилеше къ точности предпочель простоть 
самого процесса вычислевя —его самого оно не затрудняло. Между 
тьмъ примвнеше для приближеннаго интегрирован\я дифферекщаль- 
ниыхь уравнейй той формулы. которою онъ пользуется, дьлаеть самый 
процессь вычислевя гораздо болфе сложнымь а, главное, требующимъ 
гораздо болье напряженнаго вниманя, нежели по той, которую Адамсъ 
приводить, но которою не пользуется. Эта послЬдняя формула явля- 
ется осповною вь нашей статьф, и основанный на примъневи анало- 
гачной формулы премъ быль разработанъ Штёрмеромъ въ его изслЪ- 
дованяхь о сВверныхь дявяхь, хотя онь на Адамра вигдв не ссы- 
лается. Столтъ только сличить нашу схему вычисленя съ процессомъ 
описаннымъ на стр. 46—49 книги ВаяЧомь’а, чтобы видёть какъ, 
пожертвовавь незначительно въ точности, оказалось возможнымь да- 
стигнуть крайней простоты процесса. Все наше вычислеше выполнено 
цфликомъ при помощи 50 сантиметровой логарифмической линейки, — 
у Адамеа семизначными догарифмами, полученные же результаты 
отличаются оть чисель Адамса не болЬе какъ на 2 или 3 единицы 
нашего послёдняго (четвертаго посл запятой) знака, что составляегь 


менфе нежели оо миля. въ размёрф канли и значить боле чфмъ 


достаточно для разсматриваемаго физическаго вопроса. 

Веобще при производетвЪ всякихь численныхь вычислени надо 
руководствовался празиломъ;—точность вычнелен!я должна соотвфт- 
ствоваль точности данныхь и той практической потребноети, для которой 
вычиелеше производится. Этимъ правиломъ надо руководетвоваться 
и при численномъ интегрирован: уравнен!й; соблюдене его знали- 
тельно сберегаеть трудъ и время. 


$ 1. При разработьЪ конструкцЕ! ифкоторыхь судовыхъ устройетвъ 
и механизмовъ приходится имфть дЬло и съ такими вепросами ди- 
намнки и теорйг еопротивлен!я матераловъ, которые требують рёше- 
ня обыкновенныхь дифференшальныхь уравненй. 

При этомъ по большей чаети нЪть надобностл непремфино имфть 
общий интеграль уравнешя съ входящими въ него буквенно про- 
извольными постоянными, а достаточно найти частное рёшене удовлет- 
воряющее даннымь начальнымь усломямъ, и при томь для ограни- 
ченнаго промежутка измЪфняемоети аргумента. Для практики нъть 
также надобности въ аналитическомь прелетавленнг искомой функциь, 
а лоетаточно составить таблицу ея значеши или построить кривую ее 
представляющую. 
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Премы такого численнаго рышешя диффер.ур-н рфдко излага- 
ются въ куреахь, между тьмъ по сути дьла эти премы не сложиве 
вычисленя опредфленныхь интеграловт,. 

Вь настоящей замЪтьВ я им\ю въ виду паложить наиболфе удоб- 
ные и обще способы такого чиеленнаго Беня обыкновенныхь диф- 
ференщальныхь уравненй. 

$ 2. Конечно и здфеь начало положено Эйлеромъ, который уд®- 
ляеть УП главу Гго тома ШЩановез Сус ПиедтаИ$ изложен 
хпособа приближеннаго численнаго рёшеня диффер.ур-вя перваго 
порядка, къ которому всегда можно свести все дЪло. 

'Итакь положимь, что искомая функшя у опредфляется диффер. 
ур-Шемъ 
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и начальными услошями;: при х=а должно быть у=9. Выбравъ до- 
статочно малый промежутокъ №, можно принять, что пока х закля- 
чается между а и я--й, у мало отличается оть 6, п производная 
сохракяеть постоянное значен{е [(а,5), такъ что для этого промежутка, 
будеть: 
у—6--«—-9Раб-Ре—а.4, 
и ствдовательно для конца его. когла х=а--й будеты 
= А.йЬ. 


“Такимъ образомъ мы получили приближенное значен!е у=6, соотвфт- 
ствующее значеню ха, =а--й. 
Совершенно также полагая 


А, =Иа,.5} 
получимь для 2=а,—а--28 величину 


в, = --А,.й 


инт. д такниъ образомь составится такая таблица соотвытетвующихь 
значен\ 
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при чемъ 
4.— На», в) 
и 
Фа -еа, В. 


Изложивъ этоть’‘иремъ, Эйлеръ говорить: „чфмъ меньше промежутёи 
между послЪдовательными значещями х ТВмъ точиЪе опредЬляются 
ВсЁ остальныя величины; тВмь не мёнзе велЪдетвье большого числа 
накоплеше и этихь малыхь погршностей можеть достигвуть значи. 
тельной величины“. 

„Погрышность при этомъ вызислеши происходить оть того, что 
на протяжеи каждаго отдфльнаго промежутка об перемвиныя х 
и у принимаются сохракяющими свой постоянныя значеня. соотвЪт- 
етвуюцйя началу этого промежутка, такь что и значеще фунеци 
и Иа, остается постояннымъ, поэтому чёмъ быстро значение этой 
фуньщи изнфняется при переходЪ оть одного промежутка кь слфдую- 
нему, твмъ бблыпую можно ожидать погрышноеть. Это невыгодное 
обетоятельство имфеть мЪФето обыкновенно тамъ, гдЪ значеше {(а:) 
или уничтожается или же становится безконечнымь“. 

Чтобы устранить это „невыголное обстоятельство Эйлерь посту- 
паеть такь: положниъ для простоты письма, что /(жу) обращается 
въ нуль при т=0,у-ф тогда надо положить 


в=а--ь у-=6-п 


и ваше уравнеше приметь вилъ 
=Да-ь а ....... М 


и начальное услоше будеть: при Е=20 должно быть и=0, вмЪетв съ 
тВмъ при Е==0 будеть и | )= На, р=о 


поэтому положивь 


па... | 
подетавляемь величину [3] вь уравнеше [2] и опредъляемь показа- 
тели н, р... н коэффищенты а; в: ... такъ чтобы уравнене [2] 


имЪфло мфето тожественно, иными словами по извфстному прему раз- 
агаемь величину 4 вь рядъ по возрастающимь степенямъ =. 
Въ томъ елучаь, когда /(а,5) = можно уравн. [1] писать въ вядъ 


ЧЕ 


-артеру- ВЕ 
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функшя Реп) при ==0 и п=о будеть обращаться въ ноль и зна- 
чить надо будеть взять разложеше вода: 


аз рале-Е ыТ-- 
и вайдя коэффищенты В,,В,........ и показатели А, в....-.. 
этоть рядъ „обратить“, т.-е. выразить ч въ вндЪ ряда расположеннаго 


по степенямъ =. 
Такь Эйлерь береть такой примфрь: 


©) 
и при х==а должно быть. у-=а. 
Положив 
д—а-ь, у=а-рт 
позучимъ 
ЧЕ И) НШ... 9 
дФлая затЬмъ: 
деи" На +. + 
имфемъ. 
ата ти "+ раде -- я 
24] — паи На] 52. Че о. ] 
Откуда ельдуютъ уравненя: 
й—1=1 
которыя дають у . 
ыы 
и значить 
И ЗИ ИВ 8 
ы в а О Ы 


Легко видЬть, что если положить ===—й, то п будеть разлагаться въ 
рялъ вида 


п АЕ Ве Св... .. ы} 
и выфого того, чтобы находить 4, В, (,....... череть обрещене 
ряда 45) проще вернуться къ самому уравневю (**} и подетавивь въ 
него выЪето $ и п ихь величины искать д, В, С,...... что деть 


пы 


ща 
или 
ау в—9 411 ев: +. 


значить 2 должно быть меньше д, чтобы у оставалось вещественнымъ. 
и при 2=0,—@—й будеть 


: 
ув, —-а-Руай-- у ты + Г №: 

Пару =: =в— 1; у=0, можно теперь принять за исходную и при- 
мБнять вышеуказанный процессъ. 

$ 3. Чтобы устранить первый изъ указанныхь имъ недостатковъ, 
т-е, что функшя Ге,у) привимается для каждаго промежутка 38 
постоянную, тогда какь на самомь дЪшВ она претерифваеть изуфне- 
я, вь особенности если промежутки не весьма малы, Эйлерь пред- 
полагаеть, что въ смежности съ каждымь изъ разоматриваемыхъ эна- 
ченй перемфнной независимой искомая функШя у разлагается по 
Тейпорову ряду, такъ что будеть веобще: 


уе О ум. 


выражешя же произведныхь составляются пользуясь даннымъ урав- 
нешемъ 


УИ... Ч 
дифференцируя его послёдовательно п замфняя у’ его величиною, 
такъ что будегь: 

НИ 


арын абы" — Ри ЕРИР ЕЕ 
Полагая въ этихь выражейяхь 2=а,,у=0ь. найдемь зкачешя 
И, Ч" и... - и двлая д—0.:==0, | по фор. (4) найлемъ 
У=Фь+1 т.е. по парё (4,0) получимь слфдующую пару (а +: 9+1) 
Такимъ образомъ все дфло сводится къ вычиелензю количества: 


6. .—-Ш=4. ВВ й ес," В+... -® 


иными словами „кооффищентовь“ Ак, Вь,Сь (о --.--- которые 
находятся или по фори. (5) или какъ замфчаеть Эйлерь „и не прои- 
зводя дифференцированя“, а при помощи подетановки 

у=Ы-И РЕ 
вь уравнене у’ =/(т.у) п разложев я величины въ ряль по степенямъ 
Б, чему Эалерь и даеть затбиъ примЪъры. 


8 4. Изложенная метода Эйлера была затфмъ развита Коши, 
который начинаеть съ того, что по заданному дифференшальному 


уравнению 
УИ.) ©) 


и начальному услов: при х-=л, должно быть уу, составляеть 
рядъ равенствъ: 

= (и —Ихь о) 

(&—#) Иа, у) 

— 2) И, 9) 


Уи. (Х я, ОА, 1.41 
и береть ихь сумму: 


уу) арфа, - 

Ч Иа: ® 
т.е. дЬлаеть то же самое, что и Эйлерь, но затЪмъ Коши доказываете. 
что при соблюдеши нЪкоторыхъ ограничительныхь услов относи- 
тельно величины промежутка Х—х и функци Ге) разсматриваемая 
сумма, при веопредъленномь возрастании: числа я и уменьшен каж- 
дой изь разностей х,—х,_., имфеть конечный и опред®ленный пре- 
дьлъ, представляюций такую функц перемфнкой Х, которая удов- 
летворяеть уравнею 


ЕХ=ИХ.ЕХ)]....... 8 
равносильному ур. (1). КромЪ того, онъ показалъ какъ найти предьлы 
погрышности, которую мы сдфлаемъ если промежутокь Х—2 подраз- 
дЪлить иа какое нибудь конечное число частей, иначе когда разности 
т... -- будуть имфть конечныя значеня, какь это и 
сть въ методф Эйлера. 

Но вычислеше этихъ предфловь погрышности не представдяеть 
практическихь удобетвь, поэтому доказательство Котни имфеть глан- 
ныхъ образомъ теоретическое значенще, какь первое строгое доказа- 
тельство того, что всякое дифференшальное уравнене перваго порядка 
нмЪфеть общий интегралъ, 

$ 5. Пикаръ въ первомъ том своего Туайё 4 Апейузе, изложивъ, 
вышеупомянутое доказательство Санейу существовашя общаго интег- 
`фала, даеть свое собственное докалательство, основанное на примфнени 
особа нослфдовательныхь приближенН!. 

Сущность способа Пикара состоять въ стёдующемь: положниъ 
что дано уравнен® 


МЕ... :..: 0 


и начальныя услов1я: при х==х. `должяо быть у=у. 
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Вм%ето этого.уравневйя лля перваго приближеня Никарь береть. 
такое 


у" Их, 4) 
откуда слфлуеть 
а 
о 


Для второго приближен1я онъ полагаеть 


уз’ == Да, у) 


такь, что 


= (ле, 4-ю 
г 


вообще #11)" приближене получается по предыдущему, взявъ вмЪ- 
сто предложеннаго уравнешя такое 


==, д) 


« 
у Лова-ь де № 


о 


Пиваръ затЪмъ доказываеть, что когда функшя {}(2,у) удовлетворяеть 
иЪкоторымъ весьма общимъ условямь и разность ;— не превызна- 
еть НЪФкоторой границы, то описанный процеееъ есть процессь еходя- 
зися, такъ что по мЪрЪ возрастангя значка Л, функ я у. етре- 
митея кь яъкоторому предёлу Р(1), для котораго будеть 


} ИР, 


45 


Е&а) . 9 


т.е. величина, 
у=Е(9) 


удовлетворлеть предложенному уравненно и при х-=х, принимаеть 
требуемое значеше у». 

Уравнене [| заключаеть въ себф критерШ, по которому можно 
еудить при практическомь вычиеленлт достигнута дли прелфльная 
фунешя т) или ибть, и если нЪъть, то велика-ли погрыашноеть, 
такъ какъ покл предьль не достигнуть послълующее приближене 
зе будеть совпадать съ прелшествующимь’ ему. 


—_п- 


Пикаръ развиль евою методу ие вь цЪляхь практическаго вы- 
численя, а вабъ мегоду дающую строгое доказательетво существовал 
шя общаго интеграла. Эта метода одинаково примфнима и къ систем 5 
обыкновенныхь уравненш, такъ напримфрь для снетемы трехъ ‘урав- 
нений 


= (,н, 5,6) 


и начальныхь условии: при #= 2% д 
слфдовательныя приближеня находятся по формуламь. 


„я 


) но -- то 


ан 


о 
ель НЫ - п 
^) 


+ 


1 -| Диона ку 
2% 


$36. Приведенный выше процеесъ поелфдовательныхь приближе- 
НЕЙ остаетея сходящимся и въ томъ случаЪ, если принять за походиое 
не первое приближеше по Пикару, а какое либо иное ботъе близкое 
напр. получаемое по методь Эйлера и чВмь это ‘приближеше ближе 
кь истинному тёмь быстрЪе будеть достигнута предьльная функщя 
(2) съ избранною степенью точности вычиеленй. 

Но еели брать исправленную методу Эйлера, т.-е. ту, гдЪ оНъ 
береть не только первый, но и сльдуюние члены разложешя въ Тей- 
зоровомь ряду, то практическое выполнеше вычиелений часто оказа- 
лось бы длнннымь и потребовало бы много излишней работы. Чтобы 
этого избъжать было предложено два способа, изъ коихь одинь при- 
задлежить проф. Геттингеискаго Университета Ваззе, другой проф. 
Универсытета въ Хрисман УМомтему. 

$ 7. Метода Рунге въ примфнеши къ уравненр 


У ее 0) 


состоить въ развитш приведенныхь выше словь Эйлера: „такъ какъ 
Ебе ДЬло сводится къ нахожденю количества 
в д 


даб Ах. ЙЕ Вь ИУ [12} 


—_ 8 — 


экю можеть быть сдълано и не производя дифференцированя*, послЪ 
чего Эйлерь примфняеть подстановку и разложеше въ рядъ при по- 
мощи неопредфленныхь козфищевтовъ. 

Въ 8 $ приведены обийя выражешя производныхь У” я У” по- 
лучаемыя дифференцированемъ предложеннаго уравненя; полагая 
въ нихь 


&=@ь Я=ь 
А, = Па») 
Вы а, 4 Аба) 
Сы (вы) Ау") -- Зе, Вир ава 
или положивь для краткости 
(ав дь) = Рь ; Каф) = 
Тиба беты Раб) яв аб) 


получим 


и 
Пи (вы) -Роьрыи (в, бы) -- А ав 6) = Рь 
Ву’ (в, 6) 
будемъ имфты: 
В Р,-Адь 
Сь=РьЕ®ь 


Изь формулы (12} видно, что если при вычнелеши разности 

-.—в: ограничиватьел членами третьяго порядка относительно д. 

то нЪть надобности вычислять въ отдфльности величинь 4,, В,. С,. 
а достаточно вычислить сумму: 


АрЕаВИ-ЕОвЕК,Ь 
при чемъ 
а авАНтом. 
Чтобы вычислить величину этой суммы или, что то же, множителя 


№, обратимся еперва къ чертежу. 
Положимъ, что координаты точки М суть 


х=Об,=в: у=беМ-АЬ, 


к ММ, есть истинная интегральная кривая, проходящая черезъ точку М, 
т.е. кривая, ордината произвольной точки которой, разсматриваемая 
кавъ функийя абециесы 2. удовлетворяегь уравнею 


у—=Ах,у.. 


— 9— 
Вообразимъ касательную къ этой кривой въ точкЪ М и отыфтимъ 


точки пересбчешя Н, и Е этой касательной съ среднею ординатою 
Е, С н съ крайнею @,.,М, 


У 


о] % 


Фит, 1. 
Очевидно, что тангенсь угла наклонешя этой касательной есть 
4,=Ка, 5} и ординаты точекь Н, и Ё суть 


НЫЕ ТАЙ и Ч. -ЬЧАА 


Нодставляя во вторую часть уравнен:я (1) вуфето с величиву 
т и выфето у величину п» получимъь уклонь касательной къ 


интегральной кривой, проходящей черезъ И, въ этой точкф, обозна- 
чимъ тангенсь этого уклона черезъ а,, такь что будеть 


ва ь-УАА) 


преведемъ прямую Я,К подь этимь уклономъ и отыётимъ точьу ея 
пересьчешя К сь ордиватою Сь,,М,. Такъ какъ точка М, весьма 
близка къ точкВ С (разстояве Н:С второго порядка относительно #). 
то прямая Н.К будеть приблизительно параллельны касательной въ 
точкЬ С къ интегральной кривой, проходящей черезъ М, значить если 
провести черезь точку М прямую наразлельную Н.К, то точка ея 


— 80 — 


пёресфченя У будеть лежать весьма близго къ М, во всягомь спучаь 
гораздо ближе точки РЁ, соотвЪтствующей первому Эйлеровскому приб- 
лиженю. ВмфетЪ съ тьмъ очевидно, что точка С интегральной кри 
вой будегь лежать между точками Н, и р и при томь весьма близко 
кь середин® между ними. ПослЪ этого общаго объясненя качествен- 
ной стороны дЪла сославимъ выразжевя ординать точекь М, Я, ри Л 
и уклоновъ касательныхь къ интегральнымъ. кривымь, проходяжимъ 
черезъ эти точки. 
Очевидно, что координаты вашихь точекъ будуть: 


Ма, ь; Н, ( а.ть ы-т 44) 
Дань ый); (азы). 


Обозначая соотвфтетвенно черезь т, и &, тантенсы уклоновь ка. 
сательныхь въ точкахь 0 и +) будемъ имвть: 


1 1; 
= Цыуьь +) 
== (аз-- . -й). 
Разовьемъ теперь величины а, т» ®, разлатая ихь по степенямъ 
буквы й п удерживая въ этихь разаскещяхь лишь члены до 2-го 


порядка включительно. 
НуБемъ: 


К д ки 


Балви] = ь+ 3. ВА 1 РУз 


те Даь БЭН Ма оду аь Ба Е рава 


норы во о ра, БД ев 


О И И 


ивр Ва ВР. 1, 


Ё . | 
Замфнивъ въ выражещи т, величину эй черезь.й, получимъ 


БА ВАР. 0. 


—_ Я — 
Намь надо составить величину 


№ 


т „1 
ВЕР 
«топить толька взять выражен! 

г [анны] 


аналогичное тому трехчлену, который составляется изт, ордивать ири 
вычиелени площадей по формулЪф @имиеона, ‘чтобы убфдиться, что 


ы т 
м ааыы (о... 


и. значить. ордината /,.: будетъ 
Выль 


Итакъ. процессъ перехода оть пары значенй са; у==Вь въ слЪ- 
дующей парф. л=а =а.--й: уУ=А, а вынолняется слфдующимь 
образомъ: 


1) Вычисляетея величина А, по формулЪ 
Ара, сео. 04) 


2} Вычисляются величины а п 3, по формуламь 


ве день там | 


- 15) 
ыы Да А, ы-роь®) | 
3) Составляется множитель №, по форхуль 
м [анны | и. 0 
4) Величина 6.1 будеть 
мА, с... 


Это и есть метода Рунте. 
Замътимь. здвеь же, что молелно вычиелять еще н величину 


пе вы-то. 2.09 
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чтобы имфть контроль, а также ординату и уклонъ касательной въ 
промежуточной точеЪ С кривой, чтобы имфть возможность сразу вдвое 
уменьшить промежутокь при исправлеши интегральной кривой по 
методЪ Пикара, какь это можно будеть видфть на чиеленномъ при- 
мЪрь приведенномъ виже. 

Главный недостатокъ методы Рунге соетоить въ томъ, что при- 
ходитея вычислять по крайней мЪрё три сери значейи функиш 
ху) и значить, когда предложенное уравнене сколь нибудь слож- 
ное, то эти вычислешя потребують затраты большого труда, въ значи- 
тельной дол напраснаго, какь это будеть видно по прим Бру приво- 
димому ниже. 

8 & СЪ 1890 г. по почину Норвежекаго проф. Биркеланда быль 
снаряжень рядь экспедищй въ полярныя страны для производства 
енстематическихь наблюденЕ! надъ сфверными счянёями и магнитными 
возмущен ями (бурями), чтобы точно установить связь между этими 
явленями, которая была подивчена уже давно. 

Въ поелЪдинхь экспедищяхъ принималъ учасце и проф. К. Итёр- 
меръ, который разработаль математическую теор сБверныхь ся, 
вполнЪ объясняющую вс наблюдаемыя язленя. По этой теорш сЪ- 
верныя офявя образуются назлектризонанными частицами (чаотицами 
электричества) выбрасываемыми солнцемъ, въ особенности тьми извер- 
женлями, которыя намъ предетавляютея въ видь солнечныхь пятенъ. 
Эти частицы, двигаясь съ громадными скоростями, попадають въ еферу 
магнитнаго поля земли и достигая верхнихъ, весьма разрёженныхь 
слоевъ земной этмосферы, вызывають свфчеше, составляющее свверное 
чине. съ другой стороны онф производять подобно электрическому 
току магнитныя возмущеня. 

Для установлевя этой теор ИНёрмеръ предириняль громадную 
вычислительную работу — именно опредвлеше вида траекторий! опиеы- 
ваемыхь разнаго рода наэлевтризованными частицами въ магнитномъ 
пол земли. 

Траектор{и эти весьма разнообразны и сложны, предетавляя кри- 
выя не плоещя, а иногда и спирали съ нфеколькими завитками. Лиф- 
ференщельныя уравненя движеНя частиць не только не интегри- 
руютея въ конечномъ вид, но не интегрируются и въ квадратурахъ, 
и Штёриеру пришлось прибътнуть цъ приближенному способу чиелен- 
наго интегрированя. Для этонй цёли онъ развиль особый премъ, по 
цоторому ихъ самимь и его асепстентами было вычиелено около 150 
траекторй разныхъ видовь для изученя ихь типовь, при чемь для 
наждой траевтори приходиловь вычислять по 190 — 120 точекъ, т. ©. 
:0 810—369 координать. Эта работа потребовала около 5000 рабочихь 
заеовъ и, ельдовательно. метода Шиёрмера въ достаточной мЬрЪ про- 
върена на прабтнкЪ. 
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КромЪ того въ осповЪ евоей она весьма близка къ тВмъ мето- 
дамь, которые прямфпяются въ небесной механикё подъ вазвашемъ 
„методы механическихь квадратуръ“ для вычисленя возмущен пла- 
неть и комегь и которая провфрена столЪтнею практикою спещально 
вычислительныхь учрежден какъ при главныхъ обсерваторяхъ, такъ 
и при издан!и астрономическихь таблицъ п альманаховъ. 

Мы изложимь методу Штёрмера не слЪдуя буквально за, ея авто- 
ромъ, а ладимъь болфе элементарный и оби выводъ его основной 
формулы, указавъ опять таки ея связь съ методою Эйлера. 

$ 9. ПШтёрмеръ развиваеть свою методу для уравненй второго 
порядка, не содержащихь первой производной неизвфетной функфи 
и для системь такихь уравненш, мы же начнемь сь уравнений 
перваго порядка, и вмЪфсто того, чтобы нриводить, какъ то дълаегь 
Нибрмеръ, велкое уравнене къ системф уравнен! второго порядка 
разсматриваемаго пмъ вида, мы будемъ приводить ихъ кЪ систвив 
уравненйй перваго порядка. Но кромф тото мы покажемъ отдьльно 
и еамую методу Птёрмера для непосредственнаго инхегрированя 
уравнешй второго порядка какъ наиболье важныхь для вопросовъ 
дивамики. 

Итакь положимь, что по прежнему предложено уравневще 


УИ о. 


и дано начальное услове при х==а должно быть у==6. Мы видёщи 
при изложенш методы Эйлера, что вс дёло еводитея къ вычисленю 


6, две Дори о. © 
разности между значешемь В, нензвъстной у, соотвфтетвующимь 
значеню #—=а-- = и значешемь 6... соотвЪтетвующимъь йь+1 
причемь безразлично какимъь способомъ эта разность вычиелятся. 

Такь Эйлерь вычисляеть въ отдЪльности коэффищевты А» В»ь 
Оь Ль ... составляя по данному уравненю постЬдовательныя пройз- 
водныя искомой функши. 

Рунге, ограничиваясь первыми тремя членами вмфсто производ- 
ныхь В, С, второго и третьяго порядка вычиеляеть три’ зизченйя 
производной первато порядка непосредственно по данному уравненю 
но для разныхь паръ значениях н ун по нимь составляеть требу- 
емый трехчленъ. 

Шутёрмеръ не ограничивается подобно Рунге первыми тремя чде- 
нами. его метода даеть возможность брать скольве угодно членовъ 
въ разложенфи (*), пронаводя при этомь вычиелеше лишь одной серии 
значений величинь 4х непосредственно по данному уравнению а не 
вЪеколькить, какъ того бы требовала метода Рунге и совевыь не преиз- 

2 
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водя непосредетвеннаго вычислешя производныхь высшихь. поряд- 
ковъ Вь Оь Рь .... какь то дьляется въ методЪ Эйлера, а восполь- 
зовавигись для ихь нахождешя или исключешя составлешемь „раз- 
ностей“ перваго, второго, третьзго и т. д. порядковь величинь 4»й 
Этоть послВдёй процесеъ требуеть лишь вычитаы и по своей про- 
стотЬ конечно не можеть быть даже и сравниваемь съ процессомь 
непосрелственнаго вычисленя значенй А» по ланному уравненю или 
значенш В, С, Р,... по уравнешямь получаемымь изъ даннаго 
дифференцироващемъ. 

Итакь положимъ, что составляются параллельно дв таблицы 
Ги П, изъ которыхь первая содержить послёдовательныя значеня 
неизвЪстной у н ихь первыя разности. вторая послЪдовательныя зна- 
зеня величины 


Пи АЙ == (а, вый 


и ихь разноети перваго, второго третьяго и т. д. порндкойъ, причемъ 
для престоты мы ограничимся лить третьимь порядкомъ. 

Положимь, что въ этихь таблицахь мы дошли до величинь 
соотвзтствующихь значеню х-==а -=а-|- эй такь что напи’ таблицы 
будуть такого вида, 


| 
| 


} 
} 


Разности випонваютея въ промежуточную етрову между тёми 
двумя числами, для коихъ они берутся такъ 

Дивы 

Ат, 

А» 


ДВ» итд 
Ци 2— Ц. : Дт Пи Пи ит. 
Ап». — ДАць-з : Ан _з=в Ди: — Аи 


Весь вопросъ и состоитв въ томь, какимь образомь добавить къ 
этимъ лвумъ таблицамь по одной косой строе снизу т.е. перейти 
оть и вь в 
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Мы этого достигнемь, если получимь какимь бы то ну было 
‘способомъ величицу 4$, ибо прилавъ ее ть ё, будемъ имЪть 


= 5-Е Ав, 


и первая таблица будеть заполнена. Вифесть съ твмЪ сейчаеь же 
прибавится одна косая строка и ко второй табдиц®.—въ самомъ дфлЪ 
по данному уравненю 

==, 


полагая въ Нем 2—4, У Вы 
находимь 
на=- (ани 


и значить 
Аи = Канны, би) 


и внисавь пи: ВЪ Свое МБОТО сейчась же получаемъ 


Ап, пу: — 1. 
и затвиъ 
А? Ал, — Ап,.-13 Ат, = Аи, А,» 


и 


и сльдовательне вторая таблица будеть заполнена и можно будеть 
перейти къ елфдующему значеншо я--2 и т. д. 

Такимъ обраломъ надо только показать какь вычислить Ай, по 
имфющимеся въ таблицЪ Н чиеламъ, т. е. по 


пл, ть в, Ацы А а-- 


Это производится по стВдующей формулЪ: 


А+. 8) 


А. Ай» 


которую мы сейчась и выведем. 
Цо Тейлорову ряду согласно сдёланному обоаначеньо имземъ: 


ААВ не" 


или опуевая для простоты письма значекъ ю въ правой части 


ом. 
ии СТИ 


ис мы пиъемъ 


Е 
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слфдовательно, но тому же Тейлорову ряду будеть: 


= вс ОИ а 
Пи — 4 — ВА? -- в ГЕ 
вом, А 
чей —2 8-40: 807-168; 
[о 28 


в А „ 
мы 48—85 9С"— 270-81 


Но по самому соетавлетю разностей 


Айни Ин; Ап эти — Э.Н; 
Аут ==» — ЗЧя-1-Е ЗИ — Ив с. 
на основан этого будеть 


РН № 
А.В, 50, Е, 


2+.-- 


В 28 и 45 
ДА?» з= 2—6. риЕЫ— 


АЗциз= 


Намъ надо составить величину 


Г 23 р [2 
ды, -в" С НЕ, |... 
Чтобы этого достигнуть, стоить только умноживь предыдуния 
уравнешя на множители а, 8, т. прядать ихъ къ уравнению 
Пи АЙ 
и опредфлить неизвЪетные множители танъ, чтобы равенство: 


пиар, ВА ти: = Ай - в +0 2..6 г |... 


нивло мБето до членовъ возможно высокаго порядка отосптелино №. 
такимъ образомъ для опрелЪленя а, В, т имвемъ уравненя 
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которыя дають 


и значить будеть 
1 БИА ла . 
Др —ты-Н о Ато Рада Ро: ААМь-ь (12) 


Само собою понятно, что эта формула лишь приближенная и 
точна только до членовъ 4-го порядка относительно Ъ, ибо разность 


2 
120 31 


ыы = 


а не нулю и, слфдовательно, первый изь отброшенныхь въ иреды“ 
дущей формулЪ членовъ есть ЕЕ. 


Не трудно видъть, что если бы мы хотфли имфть для Ав, фор- 
мулу вЪрную до членовь 5-го порядка, то надо было бы принять и 
разноети 4-го порядка, т.-е. Аи; и если бы хотЬли имёть п члены 
6-го порядка. то надо бы принять и Дб; и т. д. и наша формула 
была бы вида 


1 5 5 
Альт АА Аи Ань р еАУнь (09 


ибо очевидно, что кооэффищенты о, В, т опредьлялись бы прежними 
уравнен!ями и сохранили бы свои величины и добавились бы члены 
<ъ коэффищентами 5 и е при чемъ оказалось бы 


Наобороть, если бы мы предпочли брать болфе мелые шрюме- 
жвутки № и довольствоваться лишь членами 3-го порядка, то при 
вычислени можно было бы ограничиться разностями второго порядка 
и наша формула была бы 


1 Эла 
АБ, эп + А 
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и вь этомь случа таблицы приняли бы такой видь 


т п. 
И ПЗ М ном 
—— - м 
#—2 | в. | и | 424, 
2, Аи-+ | да 
|: Ь т | Ау» 
А Ал, 
Юз р [о В | 


а | 


при чемь точность, ими дославляемая. одинакова съ даваемою методою 
Рунгс, простота же вычислешя будеть видна на приводимыхь ниже 
лиеленныхь примфрахъ. 

$ 40. Въ предыдущемъ локазано, какимъ образомъ яродолакеиль 
таблицы, когда верхёя строки ихъ заполнены и мы имфемъ въ ниж- 
ней косой строкЪ всЪ разности до требуемаго порядка, необходимо 
показать какимъ образомъ начать составлеше атихъь таблицъ когда 
для "=0 будуть извЪетные лишь чиела 


и непосредственно по нимъ вычисляемое 


П==й. Аа). 


Яено, что если бы мы хотбли нести вычислеше дю разностей 
3-го порядка то намъ надо бы получить числа соотвЪтетвуюнщя зна- 
чемямъь и=0, 1, 2, 8 т-е. имть для начала таблиць: 
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Этого можно достигнуть двоякимь путемъ: 
1) Положивъ въ предложенномь уравнени 


мы 


разложить по извфетнымь правиламъ 8 въ рядъ ио степенямъ фи 
положивъ затЪмъ $-=й, 27, 84,... вычислить соотвытетвуюция значеня 
у, которыя и будуть представлять 4.. 6». 6; Цонятно, что здфеь само 
вычислеше покажеть какой величивы надо брать промежутокь Л, 
чтобы отбрасываемые въ ряду члены не вщялы на результать при 
избранной степени точности. По полученнымь такимъ образомъ па- 
рамъ значений (а 5,), (а, 5.), (а, 9.) вычисляются п,, Ц» цз, ПО БИМЪ 
ихъ разности и таблица будетъ подготовлена для продолженья. 

2) ПослЪдовательнымь приближешемь пользуясь разностями не 
нижней, а верхней косой строкв. 

Мы имЪли общую формулу 


1 5, з 
петр Аа иь о-в Аа 


Считая, что третьи разности постоянаня, т.-е. 


ДЗ Ан ата АЗ, 
буденъ имЪть 
А? 


=. —Д*и,- = А» — 249 


Ал»-: =Ал„— АЯ, ‚= Ап— Аи Аа 
к значить будеть 
1 В 1 
АБ, А -—15 Рпи-- Аи (9 
1 5 до Е Ах ый 
пе дл, А, А о с 9 


Аа о В, А, 


Полагая въ фор. () я=0; въ фор. (*% я и вь фор. (°) я 
получимъ: 


то АРА А... 64 


2 
Ав орви Ам. ... 


др-ь--- Ао АыВ ,. .. 09 
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Если же ограничиться членаии второго порядка, то будеть 


1 т 
АВ ты Аж 0... 7 


А А-а Аь ..... 


Посл довательность приближешй ведемъ слЪдующимь образомъ: по 
даннымъ а, и $, вычисляемъ п» затВмъ беремъ 


4—1 
= Ат, 
по полученной парв (а, 6) вычиеляемъ п, п беремъ разность 
Ат =: — 
вновь вычисляемь ДЬ, по формулЪ 
Ап Ат, 


и новое исправленное значеше $,-—5-- 5, по этому значеню Ь, вы- 
числяемъ ц, и затЬмъ 


1 
Ав, =. Ат 
и находимь В, =6, -- 46, и по полученной пар (я, 6,) вычисляемъ п, 


найдя которое составляемъ разности и находимь Ато п А\., 
Пользуясь этими разноетяин вновь вычнеляемъь 


1 1 
Ат Аль — 5 АР 


А 
Га.) 


вы 
" 


1 5 
Аню А у 


В, -- АВ, 


Получивъ эти неправленныя величины, если хотимъ идти до разно- 
стей 3-го порядка то вычнеляемъ сперва Аб, и В, по формуламь 


1 5 д. 
А, =т,---; Аи, Ге Ато 


=, АВ, 


— я — 


и вновь перевычисливъ 6, 8: В, по фор. (14 — 16) повторяемъ это вы- 
числеве, пока два послЪдовательныя приближеня не совпадуть въ 
предфлахъ избранной точности. Здфеь’ выгодне, чФмъ идти до раз- 
ностей 3-го порядка взять сперва болфе мелье промезнутки #. такт, 
чтобы разностями 3-го порядка можно было пренебречь и.ограни- 
чиваться фор. (17—18). 

Когда по ходу вычисленя будеть видно, что разности второго 
порядка становятся настольго измфняющимися, что третьи разности 
уже оказывають вщянще на результать, 16 надо или принимать ихъ 
въ разечеть, пользуясь соотвфтетвующей формулой нли уменьшить 
промежутокъ, наобороть, когда разности второго порядка огажутёя 
‘настолько малыми, что не вщяють на результатъ, то можно увеличить 
промежутокъ. 

$ 19. Само собою понятно, что изложенная метода примфнима гЪ 
«истемв нЪоколькихь совокупныхь уравнен!! перваго порядка. 

Напр., пусть дана система 


9-е) ] 


$ (2, 1,2) | 


; тогда положивъ 


. 6“ 


и начальвыя услов!я: при 2==а должно быть у=%, 2: 


пу=пий 


зоетавляемъ не одну пару, а двЪ пары параллельныхь таблиць, какъ 
погазанныя выше, при чемъ одна пара заключасть значейя 


му, Ду: ч,Ап, дэ, 
другая: 


2, Аа; в, А.А; 4%, 


ать ур. {**} а для вы- 


причемь для вычиеленя величинъ п и Е слу: 
числешя ДЬ, и Ас, формулы: 


1 5. з 
У | 


{#) 
1 5 Е 
ВА НВ Ая, ан, | 
Совершенно также вадо поступать для системы трехъ и болышаго 
числа уравнений. 
$ 12. Какъ иавфетно, венкое обыкновенное дифференщальное урав- 
нене любого порялка и всякая ‘система такикъ уравненй, если толь 


—_- 


они разрынимы относительно запвысшей производлой въ нихь вхо- 
дящей, могуть быть приведены къ систем уравнен! перваго порядка 
вида: 


Чт 


Я 


ть... 2.0 (#=1.2,3.....№) 


Такъ, напр.. если бы вопросъь шелъ о движент точки въ плоскости 
оу, то мы имфли бы систему: 


2. ас 4у 
ЛЕ Иа : у) 


4& ау 
ори, 


Вводя неизвфетныя фунюши вн х приводимь эту систему къ виду: 


#5 Чу _ 
= чи = 
Чи р 4 


ЧЕИе ий: пе о уд 


кь которой уже пепосредственно приложима изложенная выше 
‘метода. 

$ 33. Штёрмерь развилъ свою метолу въ примфненш къ систем, 
‘уравненй второго порядка, вида: 


= Ру) 


9, 4,2,0 | (А) 


ат 4,4, 5,6 


относящейся къ движению малеральной точки въ средЪ не оказыва- 
вицей сопротивлен1я. 
Разсмотримь сперва одно уравнене вида 


не содержащее первой производной неизвъетной функши, х, и воло- 
жинъ, что это уравнеше требуется интегрировать при начлальныхь 
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. о, 
условяхьы при #=% должно быть 2=ж и = =. Положивъ 


йа’ составляемъ таблицы такого вида 


канчиватьея косыми строками, водержа- 


Эти таблицы будуть 
щими чиела 


Аль, Ата в 
Ара А 


Чтобы имфть возможность продолжать эти таблицы, пренебрегая раз- 


ностями 6-го порядка, т.-е. прибавлять къ нимЪ по одной косой строкЪ 
Штёрмеръ примъняеть формулу. 


Г. 1 . 
АУ [ара Заь- |. . . 9) 


зычисливь по которой Дх,., вписывяемъ его въ свое мЪето, затЬмъ 
имБемъ: 
Аль Аа, Е Аи 
ди -- Ах; 


иелучивь величину 2,.. вычлеляемь по данному уравненю, едълавь 
въ немь #—1, х=1„., величину 


И” ЕЗИВ 


и.вписавь ее въ свое м5сто. беремъ разиоети и такимъ образамь нии 
и косыя строки таблици будуть занолнены и, пользуясь фор. (20} 
можно продолжать процессъ дальше. 
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Формула (20) выводится совершенно подобно тому, какъ фор. (12). 
охраняя прежнее обозначене производныхь неизвфетной функщи, 
въ данномъь случав х по перемфнной независимой {, и принимал 38 
‘исходное значене х„ соотвЪт. #=&-рий получимъ по Тейлорову ряду 
ДЛЯ 2. и 2%. разложен!я: 


Г ГЫ ГЫ № Га 
аа А-В С-5 т АТ 
Г [а Г > Г 
+ м 
2.10, Ай -|- В б%-+Р>, ЕР Ртэу .. 
Визств съ ТЬМЫ и, В. 
Отсюда слъдуеть: 
а м в 
ДЗ ==: — 22-9, 1 НОЕ ... @ 


Напиеавъ затЬмъ разложеня величинь 2.3.1.25 @ именно: 


г 15 #8 
ыы ЕО, Е Ри 


, м Го 
= эР--+Р- Е +18. 


и т. д. 
хоставляемь разности 


А ыы аь ар ВЕ аа 
‘= тв 5 6 24 
‚у 6 

о Е И Е ВЕ Е Ты 
те 
А ты бе Ч Расе а 


неключаемь (проще всего по метод’ множителей) изъ ур. (*) и (*) 
величины Р, Е, Р и получаемь формулу: 


. 1 Й 19 | 
о [ера ба |. 69 


которую для удобства вычислешя лучше писать въ вид 


1 


Пь[ НА оАь,— 


54 .. @- 


Чтобы по начальнымь услошямъ начать вычислене можно примфнять 
также два способ: — разложеня въ рядь. который не требуеть по 


—%— 
яенешй, и послъдовательныхь приближен! для котораго подобно- 
предыдущему изъ фор. (20) составимъ слЪдующя: 


1. И 
Ав А-а А’ д Ав 


В 1 > 1 дз 1 1 1 
Аль А-З А тт в — 4 А‘ 
1 
даа, АМ ом 


п ведемъ по нимъ вычислене какь показано въ 8 10, т.-е. сперва беремт- 
одинь членъ, съ елвдующимь приближешемъ два и т. д. 

Понятно само собою, что когда дано не одно уравнеше, а система 
уравненй вида (4) не содержащихь первыхъ производныхь неиз- 
вЪетныхь фуньщИ, то метода примфняется совершенно подобно тому, 
какъ для уравнен! перваго порядка, при чемъ составляется столько 
паръ таблиць, сколько неизвЪстныхь фуньш! имЪется. 

$ 14. Если брать промежутки настолько малыми, что разности 
3-то порядка уже не вяють на результать, то основизя формула 
будеть 


ть... 08 


_ ваты, у) 
и таблицы будуть: 


-Т: 
Е 
| 


и процеесь вычисленя своднтея къ весьма простому. 

815. Если въ методЪ Шиёрмера слфдить за тЬмъ, чтобы про- 
межутокъ № (мБняя когда нужно его величину} оставался настолько 
малымь, чтобы послфдняя изъ удержанныхь разностей могла быть 
прихимаема за постоянную, съ то» точностью, съ которою вычиелене 
производится, то результать получится сразу съ требуемою степенью 
точности и поправокь требоватьея не будеть. 
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Но если бы мы примВняли другую методу Эйлера иди `Рунге и 
желапи бы полученный результать исправить, то, какъь уже сказано, 
его можно было бы принять за исхощное приближене въ методу 
Пикара. 

Для простоты разсуждени прелположимъ, что намъ предлознено 
одно уравнение перваго порядка 


Па) ое, 9) 


и начальное услове при х==а, должно быть у-ёб и что для этого 
уравнен!я получено приближенное рьшеше 


у=Р). 


Если бы эта функщя была точнымь рышешемь уравнешя (0, то 
имфло бы мото равенство: 


ЕД. ] 


изъ которато интегрируя получили бы: 
= 
Ри) -| И.Р... © 


Еели же значеня первой части не будуть равны значенямъ той 
функщи РР), которая входить въ составъ подъинтегральной, то это 
будегь показывать, что принятое приближене недостаточно и надо 
положить 


ви | =, вме ... 89 
и принявь за новое приближеще 
у=Рт), 


повторить этоть процессь. который и продолжать пока съ принятою при 
вычиелещи степенью точности не будеть достигнуто равенство: 


Ед | [= Ека)... 2) 


Зльеь необходимо имёть въ виду елфдующее обстоятельство: 
фунвьши 2.02), Р,(®)... будуть намъ иавъетны не аналитическими свойми 
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выражещями, а будуть задаваться или опредълаться лишь въ таб- 
личной формЪ. п значить интегралы въ правыхь частяхь равенствь 
&) придется вычислять по приближеннымь формуламъ. Ясно, что эти 
формулы должны обладать тою степенью точности, которая соотвЪт- 
втвуеть принятой въ исполняемомъ вычиеленш. иначе вычиелене 
интеграловъ будеть вносить свои погрышности въ значешя функц 
Ех), Е(т)..`и мы вмфото приближея можемъ лате удаляться 
оть истивныхъ значей искомой функц. 

Премы такого вычислешя опредфленныхь интеграловь разра- 
ботаны съ большою. обстоятельноетью въ приложени къ вычислению 
планеткыхь возмущен; нфкоторыя изъ отвосящихея сюда формуль 
приведены въ монхъ лекщяхь „о приближенвыхъ вычиеленяхь“. 
Подробно же ихь можно найти въ сочиненш Т\. Оррошег, ГебтЪиев 
таг Вафифезатитате 4ег Котеёеп иоё Рапеев. 1. И. Къ сожалЬню 
это сочинене составляеть теперь библюграфическую рЪдкость 

Я приведу безъ доказательства одну изъ такихь формуль, именно 
Лапласа, достаточную для нашей цфли, 


тей г, 
| [лево фри 


1 
з.] 8% — 4] 


Ё 19 
р 


з , ]_ 863 ля ] 
— 36 [мена] воз А “5 4% |, 
которую напишемъ для кратгости такъ 
—=Т.-Рь 


а-ий 
у. Аз, 
8 


тд 


т 3 
[ни-ни эры 


ость главный члень формулы представляющ?й результать вычнеленя 
нашего интеграла по правилу „трапещи“, Р» есть „поправка“ этого 
результата. 

Понятно, что выкодиЪе и самое вычислеше разбивать на двЪ 
части, т.-6. отдьльно по извЪфетной схемЪ {„еуммы попарно“, „суммы 
сверху“) вычислять Ть и отдЪльно Р». 

8 16. Для пояснешя изложенныхь выше способовъь продёлаемь 
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нЪкоторые примвры, выбразъ ихьтакъ, чтобы самое вычислен!е частныхь 
значений функши Ду), стоящей въ правой части уравненя 


у=Ат) 


было возможно проще ибо эта часть работы поясненй не требуетъ. 
Кавъ первый примфрь, возьмемь уравнеше 


У=у, 


и начальное услове: при х==0 должно быть у-=1 и составимъ таб- 
лицу значени у для х черезъ 0,1.. 
Эта таблица будеть, если примВнять методу Эйлера, олЪдующая: 


ТАБЛИЦА 1 


в ео Ги м 
| ЕЕ 
Г 0,0 1,0000 | 0,1000 
1 0,1 1,1000 | 0,1100 
2 0,2 1,2100 | 01210 
3 0,3 1,3810 | 0,1831 
4 0,4 14641 | ол 
5} 0,5 1,6105 | Оле 
61 06 5716 | оду 
7 07 1,9488 | 0,1949 
в | 03 2.1137 | бра 
в, | 0,9 | 0.2358 
10 1,0 | 
| | 


Процессь составлещя этой таблицы ясенъ изъ формулы: 


уу»: 


= 


Примемъ полученныя значеня за пеходное приближен, т-е. 
зьлагая 


у=Р®=у 
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вычиелимь слфдующее приближене 
д 
кую | Коха 
о 
при чемь значенм функцы #№(2) булемь брать изъ предыдущей 


таблицы. 
Вычислене расположится елЪлующимь образомъ: 


ТАБЛИЦА И. 


и т. р : т У 
- И ии 

в =— |Суммыт Суммы И о ами ГВ, | Е 
Е (х). [готарио. сверху. 1 : ЧРИ 


0,0000 | 0,0000 | 0,0160 


| 
0,00101 9,0000 ° 1,0000 


ИИ 1 2,1000 . : 
110.1 1,1000 3.1000 ° 01550. 0 и] 1! 11049 
: 23100 | | 
2'0.2 1.2100 4,4100: 0,2205} 191 | {12208 
2.5410; | | | 
1,3310, 1 6,9510 | 0,3476 | 188° 13 1.3473 
2.7951 | } , | ! 
1.4641 | } 9,7461: 0,4978 | 146 15 2 14869 
| 3.0146 | : | 
| : ат 161 10 | 56105 
| 
р оо пт 18 6 1,3095 
| 9,9962 | 195 19 $ 1.9954 
4,0925 : | 
Зулаат, | | 24 р 10; 2,1995 
9:0,9; 2.3581 | 236 | № ао 
| 
1 


10: 1.0! 2,5939 'З8,4495 | 1,6735 + ми 
5: | : 


| 
Сличая величину №2) н Ра) видимъ, что равенство 
Ва) 


още не достигнуто въ предьлахь точности вычиеленёз, поэтому при- 
нявъ у), а зватять п уУ=Ее) составляемь выражешще 


Ее - Кое 
9 
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и продолжаемь этоть процессъ, пока два поелфдовательныхь ирибли- 

жешя не будуть совпадать въ предфлахь точности вычислешя. 
Получатея результаты, приведенные вь сафдующей таблидъ. 
Такъ какъ рышене нашего уравненя есть 


хо вь поелфднемьъ столбиф таблицы приведены и эти значеня, чтобы 
ихвть сличеще приближенных результатовъ съ точнымъ. 


ТАБЛИЦА Ш 


Не а С о 

0 | 0.0 | 1.0000 1.9000 1.0000, 1.0000! 1.0000 | 1.0090 1.00000 
1 0.1 1.1000 . 1.1049 | 1.1051 1.1051. 1. 1051. 1.1051 [1.10517 
В 0.5 | 1.2100 | 1.9208 | 1.22218 [ 1. 2214, 19814, 1924 1.22130 
3 | 03 1.8810 | 1.8473 | 1.3496 лалаов | сию 1.8499 [1.34986 
$ | 04 1.4641 | 1.4810 | 1.4912 | 1.4918 | 1.4921 р 1.4918 1.49188 
5 9.5 | 1.6105 р 1.6408 | 1.6475 | 1.8488 | ин 1.6488 |1.64872 
5 | 9.6 | 1.1716 . 1.8096 | 1.8299 р 1. 8219. | 1.8222 | 1.8222 ‚182212 
т | от | 19488 ‚1.9955 | 2.0099 1 2.0187 : 2.0138: 2.0138 [2.01875 
8 08 2.1487 | 3.1099 | 2.2195 | 2 зао | 2.2256 |2. 

9! 09 2.3581 | 2.4248 | 2.4501 | 2.4582 ' 2.459? [245960 
10 1.0 | 2.5989 авт лоне 2.7158 | 3 21188 [2.71838 


Сличая 22) и Ру) видимъ, что въ предфлахъ точности вычисленя 
можно считать, что предъльная функшя достигиута и дъйствительно 
Р.(®) отличается оть истинной функщи ег не боле, какъ въ 4-мь знань 
паст запятой, за который вообще при нашемъ вычислении ручатьея 
нельзя, 


отъ примфръ показываеть, что при избранной величин» про- 
межутка первое Эйлеровекое приближене, если его принять за нсхол- 
ное для процесса Никара, потребовало бы повторешя интегрирования 
нять разъ, для получешя той ограниченной степени точности, съ кото- 
рем наше вычнелеше производилоеь. Ясно, что было бы выгодн\%е пли 
уменьлитть промекутокь й пли примьнить другой ныемь для полу- 
ченёя болфе тотнага псхелнаго пополиженя, Нонятно, что по количе- 


19 — 


ству работы уменышеше промежутка влвое равносильно съ вычиеле- 
немь по упрощенной схем Руни, т.е. для уравнены 


{ 


У= Иа. 


по формуламъ 
а, =, у») 


УК пууруай 


ыы 2 
Во, По) 
аи, 


ибо въ обонхь процессахь потребуется вычиелене одикаковаго чнеда, 
частныхь значений функц /(®‚у) а обыкновенно это-то вычнелене и 
требуеть наибольшей затраты работы. 

Вычислене для нашего примфра по предызущимь формуламъ 
исполнено въ табл. Г\ 


ТАБЛИЦА №. 


1 В 
Зомая Карл | 


Вы 


. | | : 
9 10.0 1.0000 ° 1.0000 : 0.0506 } й | 
| | } 1.0500 {1.0500 ' 0.1050 


1.1050 | 1.1050 { | 
: $} 1.1602 : 1.1602: 0.1160 
00611 : : : 


0.2 {1.2710 | 12210 


1.2891 | 1.2821 0.1282 
8 |035 1.3192 ' 1.8192 : 00675 } 
В : | 1.4167 1367 

# | 04 | 14909, 1.4909 0.0745 | 


Эт 


: | ‚ 155654 0.1565 
5 | 0.5 ; 16474 164741 0.0824 } { 
| р 173298 0.1730 
6 10.6 | 18204 : 0.0919 
й Тс ГЕН Ода 
7 101 120115: | 0.1006 : 


: | 

: В 2.2Е 1 0.2112 

8 | 0.8 | 2.2227; 2.2221 | ОЛ, В 
2.3338 {| 0.2834: 

9 09 

| 


10 | 100; 27149 


2.4561 | 0.1298 


| 5.8739 | 0.2579 


211409 } 


| 
| 
| 
| 
| 
В 


По этымъ чиеламъ видно, что полученный результать пе своей точ- 
кости соотрЪтетвуегь приблизительно Ё;(2) и значить потребоваль бы 
3” 
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для получешя окончательного еше интегрированя два раза. Такимь 
образомь исходя отъ перваго Эйлеровскаго приближен!я намъ пришлось 
бы вычислить вн значении функщи (а, у), злЪеь же ихъ всего назло 4н. 

Обратихся теперь къ методЪ Штёрмера, получимъ такую таблицу; 


ТАБЛИЦА У. 


р 
© | 0.0 | 1.000 0.1000 
| | 0.1051 0.9105 
г | 01| 11051 0.1105 0.0011 
| : 0.1160 116 
2 | 0.2 | 1.2211 | 0.1221 13 
| 0.1284 129 
з |031 13495; 0.1350 в 
0.1420 142 
4 |041! 0.1492 | аа 
- 0.1569 156! 
5 |0.5 1.6484 0.1648 18 
0.1732 ое 
$ Г0,6 1.8216 | 0.1822 | 17 
| 0.1916 191 
71011 2.0132 21 
| ; 6.2114 | 22 
8 | 08| 2.2246 } 23 
] › 0.2850 235 
‚0.3 | 2.4596 | _ . : 
| 0.2587 | 
10 109 | 2.1188 
Кромь того для начала по схемамь $ 10 имъемь 
КЪ ТАБЛИЦ \. 
Ы у | у м | \ 
, , В и 
- я т 
о, 1.0000 | 10, 9.1000 | р 
| .1000 | 0,0100 | 
1| 1.1000 | 6.1100 | 
‚ 0.1300, | 
2| 1.2100 
Г | - 
0! 1.0000 Оле - 9-1009 о 1 Ду, =0.1000-—0.1050— 
6. 55 | 
| | э. о | — 0.0001 = 0.1049 
1: 1.1049 0.1158 0.1105 $ н211 0.1100--0.0150 
10.1221 | :-9.0005 =0,1158 
_ | Ду, =0.1600--0.0052— 
. 0.1051 р ^^ 6.9004 0.1081 
0.1160 : } 4.1009--0.0052-— 


—- 0.0008 — 0.1160 ^ 
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При сличеныя съ таб. Ш видно, что приблизжене, получаемое по ме- 
тодь ИШзёруера, есть окопчательное, о чемъ можно было заключить и 


и З Аз 
а рмом по величин вторыхь разностей, ибо ярно, что член вАйь 


1 : 
ке превышаетъ 5 единицы 4-го знака и значить не оказываеть вл1- 


2 
янтя на результатъ. 'Такимь образомъ здьеь потребовалось всего вы- 
числить я значений функцш /(%,у) т.е, приблизительно въ 4—6 разь 
меньше работы, нежели по предыдущимъ способамъ. Въ этомь и со- 
етоитъ главное достоинство методы Штёрмера. 

$ 17. Чтобы пояснить еще методы Рунге и Пиёрмера, возьмемь 
примфръ, приводимый въ сочинены аббата \вихпе, Бесоля Че саки 
ЧЕ етенЫс! еб Фе сайсий пубота! составленномъ по ленщямь Коши. 

Этоть примВръ слЪдующи: 


у—Уз-Нуу 
начальное услов!е ири х=0 должно быть 
таблицу аначенй функши у оть х=0 ло х 
По методь Коши въ указанномь сочинениг составлена таблица, 
значен, между которыми заключается величнна у. 
Эта таблица елфдующая. 


ТАБЛИЦА \1. 


6 требуется составить 


0.00000 : 9.00000 
0.02108 | 0.04143 
| 0.07413 | 6.17368 
0.15127 
9.14287 
9.36624 
0.50089 | 0.61275 
0.65296 | олвитв 
0.90961 | 0.96667 


0,9175 | 1.16688 


1.18879 { 1.37687 
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Въ разематриваемомь промежутьВ имзеть мЪсто неравенство 


Ка <у<Е® 
Какь видно предьлы еше весьма широве. 

Для этого уравненя при 2==0 какъ разъ имфеть мЪето го невыгодное 
обстоятелъетво, на которое обращаль виимаше еще Эйлерь —производ- 
ная у также обращается въ ноль, производныя выстихь порялковъ 
обратцаются въ безконечность, поэтому для значе близкихь къ х=0 
изложенныя выше методы непримфнимы, и проще всего разложить у 
зъ рядь по степенямь перемфнной х (эти степени не будуть цфлыми} 
и вычисливь по этому ряду значене у для какого-либо значения г 
отличнаго оть нуля, для полученя дальныйнихь значенш примфнять 
изложенные выше способы. Чтобы проще выполнить это разломен!: 
положимъ сперва #==ё, тогда будеть 


Чу 5 = 
Чав-авуу 


Делая затьмь узеё В --т) 


и 0бозначая производную по { значкомь имфемь 


А Ра 


паев... 


в/в = та-- зв 9-10 -[-....-.] 
и для опредбленя коэффищентовъ «Вт... . получимь 


Положивъ 


имъемь 


[реалов и. ] 2 зв [2 ор ВР... 


Откуда сявдують уравнешя: 
„_82 
= 


11298 = 16а 
6182 12ба8 = 168 
13381-14095 — 167 


значить будеть: 


#5 
уе 


Такимь обрааомь нолучимъ 
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Шри малыхь значешяхь х (около 0,1) полученныхь членовь доста- 
точно, чтобы пайти у съ точностью до 5-го десятичнаго знага, вЪ 
чемъ убфдимся, вычисляя у’какъ по предложенному уравлень 


'=Уж-Уи 


такъ но ряду 
й 3 

ий 21,2 _9 

у—= +\ в 


что видно но слёдующей таблиц: 


ТАБЛИЦА \1. 


0. 0 | о о 
от 9.030893 0.49201 ' 0.49202 
0.2 0.093677 


0. 


31 
9.8 9.13020 0.07222 0.97228 
0.4 0.23788 1.16857 1.16860 


0,5 0.41337 ; 1.85004 1.35004 
0.6 9.55706 1.52096 1.52108 
07 0.717381 1.68364 | 1.68386 
0.8 0.89359 1.33973 | 1.58991 


9.5 1.08519 1.99040. 1.99067 


1.00 1.29158 2.13648 2.13078 
| 


Слвдовательно. если производить вычислен!е на пять знаковъ, то можно 
принять, что у должно удовлетворять уравненно: 

УЕ 
и ч10 при х=01 должно быть у-=0.030898, соотвЪтетвеннь чему 
'—0.4920}, 

Для значенй х менышихь 041 величниа у опредфляетея приве- 
деннымь выше рядомъ, для значенИЕ же больших 0.1 мы причфнымъ 
методы Рунге н методу [Итёрмера, чтобы имфть наглядное сравнене 
между этнмн методами. 
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Для примбнешя методы Рунге съ т8мь добавлешемь къ ней, 
которое указано въ 8 7, позраздЬляемъ иромежутонъ между я =©1н 
1-=1.0 на два участка: оть д-=0.1 до =0.3 и оть 2==0.3 10 Е=1. 


1 
На первомъ участкВ беремь =) — 0.025, на второмь й==0.05, 


тогда вычислеше расположится, какъ поназываеть сл5дующая тайляца: 
ТАВЛИЦА УШ. 


(Таблица на отдВльномъ листЬ). 


Для ировфрки принимаемь найденныя такимь образомь величины 
у за цервое приближеше Ё, (2), составляемь по нимъ функщю 


уфе уз-РУу Р.Р. 


и по нему по формуть 


Ф-[ певдм-ь 
Н 


при 4:=0.1, 5-=0.030898 находимь второе приближеее, какъ это по- 
хазано въ таблаць {ГХ 


ТАБЛИЦА 1%. 
{Таблица на отдфльномь листв). 


Дия сравнешя произвелемъ тоже вычислене по методъ 1Штёрмера 
цолучияь слБдующую таблицу (Х). 


ТАБЛИЦА Х. 


(Таблица на отдфлЛЬьномь листф). 


Сличее полученныхь результатовъ ноказываегь, что метода 
Штёрмера при затратЪ работы въ пять или шесть разъ меньшей дасть 
результать одниаковый съ иредыдущимь. 

Отеюда становится понятнымь, почему Штёрмеръ предпочель раз- 
вить новую методу, нежели примфнять методу Рунге въ общирной 
работЪ имъ предпринятой. 

$ 18. Призедемъ въ заключен!е этой статьи одно замбчаше 
историческаго характера. 

ЗЪ 1826 году вышло въ двухъ бодьшихъ томахь |1 3 сочинеше: 
А, Гожевае, ТгаНв 4е$ опеНопз еШрИдиез ег Че; пиедгаез ащеневпез 
въ которомь этоть знаменитый мамематикь даеть сводьу свонхь 
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сорокальтнихь изыекан по теори Эллиитичеекихь и Эйлеровыхь 
интеграловь. 

Во 2-очъ том$ содержатся таблицы значе! интеграловь перваго 
и второго рода, те. функщи 


Въ этихь таблицахь значентя функции: Ки Е даны лля веъхъ эначе- 
в амилитулы ф черезь 18 оть 02 до 90% при вобхь значешяхь мо- 
дульнаго угла 8 (опредфляемаго равенствомь З16==е) также черезъ 
13 т-е. онф заключають 2.90.90 = 16200 результатовъ, 

Веб эти числа дацы еь 10 знаками дня @ оть 0° до 40° и съ 
$ знаками оть 45% до 90°. Самое-же вычислене производилось еъ 1+ 
п 12 знаками. 

КромЪ этихъ заблиць Гесденёте даеть п ряль другнхь весьма 
общирныхь таблиць, вычисленныхь танже съ 12 п 14 анавами при 
чемъ вов вычислен!я произведены н\гъ самим. 

Отсюда ясно, какимъ огромныхь навыкомь и какою огромною 
опытностья обладалъ Гедепате въ производетв$ чиеленныхъ вычислений. 

Объясняя способъ составлеыя свонхь таблиць онъ развиваеть 
методу разноетей дия вычислеши интеграловь простыхь и повтор- 
ныхь. Межлу прочимь, онь разематриваеть и двукратный интегралъ, 
т-е. нахождеше функщи у, опредБляемой уравнешемъ. 


'=/® 


и начальными значенямы у) и У, 
Формула, которую онъ для этого даеть, напишетея при вашихь 
обозначещяхь такъ: 


т 1 


В [ #1). 
Ааа, На Ан — зо ово в ‹ - 


ть 
= /(т,) 


и промежутокь Й. онъ рекоменлуеть брать настолько малым, чтобы 
въ этой формулЬ достаточно было удерживать только первые три члена. 

Чтобы пользоваться этою формулою, онъ составляеть двъ парал- 
лельныхь таблицы, соверщенно подобныхъ показаннымь въ 8 13, изъ 
ноихъ одна содержить величину у и ея первыя и вторыя разноети, 
вторая содержить „вспомогательную“, какъ Гесепге ее называеть 
величину п и ея разности до 4-го порядка. Преднолатая, что въ пер- 
вой части таблицы дошли до значешя у, соотвфтетв. х==х,„ онь вы- 
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чиеляеть для второй части величины п,-: Из,-о. Пользуяськоторыми 
находить нужвыя ему разности А; и Ач,_, стояпия въ той-ие 
прямой строкЪ какъ и п„, но нимъ онъ составляеть по фор. (*) вели- 
чину Азу,-,, затвмъ: 


Ду, у, Аы-ь и уна-=у,--А9, 


Изложивь эту методу онъ говорить; „Было бы желательно подобными 
же способами и при помощи столь-же быстро сходящихея рядовь 
имфть возможность вычислять послфдовательныя значешя функции 
опредфляемой дифференщельнымь уравнешемъ 


аи 
5- Раз) 


или уравнемемьъ болфе высокаго порядка. Эта задача того же рода 
какъ и задачи относянйяея кь вычиеленю проетыхь и повторныхь 
интеграловъ, мо ея уъшеше предетавляетв 10раздо большия трудности в 
90 сисё порё мы не видим дру способа ка ею достижению, какз при- 
мнеше формулы Тейлора“, иослЪ чего онъ издагаеть: методу Эйлера. 
Но стоило бы Лежандру выразить величины А, и Ам, 
черезь имфюпцяся въ его таблиц разности, етояшая въ послёдней 
хосой строкё идушей оть п, наисносокъ вверхъ и замфтить, что отбра- 
сыване вефхь разностей, конхь порядокъ выше 4-го, равносильно пред- 
положен!ю, что разности 4-го порядка постоянны, онъ получиль бы; 


Ач, := А, = А, 
,-з-- 44, -з= Аи, _,-- А, -з-Г А, = 
— Ал, «НА, НА, 


и во подетановкЪ его формула приняла бы виль 


Г. Ё 
= [мы а Нам ам, 


Это и есть формула Шаёрмера, которая вакъ разъ рьшаеть намфчен- 
ную Лежандромь задачу именно такъ какь Лежандрь имфль это въ 
виду. 

Это замъчаше можеть служить примфромъ того какъ рЬшене, 
которое представляется, послВ того какь оно дано, весьма простымь 
п очевиднымъ, ускользало оть внимая даже такнхь математиковъ, 
какь Лежандрь. 

$ 19. Въ вопросахъ инженерного дБла ветрёчаются уравнешя, 
вида 

ув.) 
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{изгибъ балки на упругомь оеноваши). Для уравнени такого вида 
можно воспользоваться тЬмъ премомт, который Лежанлръ примфняеть 
для вычиененя четырехкратнаго интеграла, т.-е. исходить изъ формулы: 


А Ща 
а р СЛ 
Ограничиваяеь разностями 4-го порядка можемъ положить 


Ат, 5 


= Ат, 1 
и наша формула будеть 


1 
А-а А, — 


При достаточно маломъ промежуткВ й послЪдий! членъ можеть быть 
В ланъ также настолько малымъ, что онъ не будеть оказывать вляшя 
на результать, и мы ‘получимь фориулу 


Ач, -з 


1 
м в 


вычиелен!е по которой настолько очевидно, что дальнфйшихъ поясне- 
ШИ не требуеть. 

Замвтимъ еще, что при чиеленномъ рышениг дяфференшальныхь 
уравнений по методь ТШтёрмера въ томъ видЪ, какъ онё здЪеь изло- 
жена, не приходится дифференцировать функции /(7,У), стоящей въ 
правой части уравнен!я, а приходится имфть дфло лишь еь частными 
значейями этой функщи, елфдовалельно эта функшя можеть заклю- 
чать п тафе элемевты, которые задаются графически что особенво 
важно для вопросовъ прикладной математики. 


ПРИБАВЛЕ ИТ Е. 


$1 Проф. П. П. Лазарев еообщиль миф, что для редактируемаго 
имъ журнала было бы желательпо имфть еще какой-либо примВруь 
приближеннаго интегрировашя уравнешя встрёчающегося въ физниЪ. 

Возьмемь какь такой примЪръ уравнен1с, которымь опредляется 
чформа мериданальнаго сЪченя капли тяжелой жидкости зежащей на 
горизонтальной плоскоети. 

Общее уравнешю граничащей капиллярной поверхности есть, 


и 


ГДЬ 2 есть вертикальная ордината точки поверхности, р; и р, главные 
ратусы кривизны въ этой точкЪ, а* вЪкоторая постоявная для данной 
жидкости. 

Въ нашемь случа граничащая поверхность есть поверхность 
вращеня около вертикальной оси т.-е. оси 2 направленной внизъ. 
Для каждой точки такой поверхности одинъ изъ главныхъ радусовъ. 
кривизны есть радусъ кривизны мериана въ этой точьЪ, примемъ 
его за ру; второй главный радзусъ кривизны р равенъ отрЪзку нор- 
малл въ разсматриваемой точкЪ, заключенному между этою точкою 
и осью вращеня. Обозначая черезъ ’ радусъ цараллели, проходящей 
черезъ разсматриваемую точку и черезъ у, уголь между нормалью и 
осью 2, будемъ имьть 


х 
эф 


и етфдовательно. наше уравнеше, будеть 


[1 за 
я $ ии с 9 


иначе: 


98 5 
Яга. сояф 
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гдЪ 4 есть дифференшаль дуги мхерижана п 
„_# 
по 


92=р:. Зтф.@ф | в 
` } 


Фр ооно. о | 


Принявъ ф за перемфнную независимую и обозначая произведныя 
значками, мы видимъ, что нашь вопроеъ приведень въ интегриро- 
вавю слъдующей систелы уравнен: 

==. Зи | 


р: . 605Ф . . ИИ 


тДЪ 


Нетрудно видфть, что эти уравневя однородны относительно ве- 
АИЧИНЪ 2, р. Го, п ихь моязно написать въ таномь видЪ: 


(9-8) 


сльдавательно 3, ги р. межно выражать въ доляхь и, иными сле 
вамп мояно безъ ущерба общности полагать в-=1 и разематривать 
систему? 


2 =р,. ЗФ 


Хр. 615$ 


Вь верхней точкЪ капли, которая очевидно лежить на оси вра- 
щен, будеть ф=0, значить начальныя условя будуть: 
при Ф=0 должно быть: 7==0 и 2= гдЪ & ныкоторая постоянвая. 
Тасъ кавъ при беззонечно маломъь ф, будеть и х безконечно 


малое то пред (*,*) пред (*) = пред (=) =( На сльдовательяь, 
и 


ый 


при Ф-=0 будетъ: 


в ^ 
Для чиеленнаго примфра возьмемъ тоть случай когда == 1, и 
Боложавь 


1==2.й; Хи. й 


тдЬ й есть избранный промежутокь между смежными значенями ф, 
который принимаемь за 5°—0,087266, располагаемъ наше вычиелеше 
какъ приведено ниже, выполняя его логаривмической линейкой 
{50 сант. длины). 

Для получена исходныхь значенш 2 и ^ проще всего разложить 
2 и го стеленямъ ф. что выполиимъ такъ: легко видЪть что р: есть 
четная функщя угла ф, поэтому полагаемъ: 


р= 1-99? В®*--...- 


тогда, 


рзто=Ф-- (- +) 7-5 «| 


ревчф —1-- («— 5) ты 


Слфдовательно, въ внду условШ что при Ф-—0 должно быть 
п 7—0, на основан! ур. (5) ныъемь: 


2 39 (т- в) ($ а в)". 


"9+ ( 


Затьмь пмфемь 


м 


2та— маф)р,==Ф- (. 


Уравнивая эту поельднюю величпну величинь 7, получаемь для 
опредьлены неизвфетныхь козффициеятовь а н В уравненя: 


ЗЕ 
5“ зв 


т.т 
в — 
18 (24 
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которыя дають 


Подставляя, нифемъ: 


. 6) 


Подагая въ этихь выражешяхь послЪдовалельно ф-=0°, 50, 10°, 
15°. т.е. дълая ф равнымь 0; 0,087266; 0,17455; 0,26180, получаемь 
слвдующую таблицу исходныхь значенй. 


5 | 1,0038: 0.0870 | 0,9943 
10° | 1,0150 ' 0.1728 99.9780 
р 


025461 0,9527 


пыфя которыя и продолжаемъ вычислене по общей схемъ, приведен- 
ной ниже на отдЪльномЪ лиеть. 

По найденнымь велитинамъ коордпнать 2 и / непосредетвевно 
вычерчивается. вынеся ва графикъ соотвфтствуюния точки, форма 
меридтана; если же кромф того желательно имЪть и его яволюту, то 
предиочтительное не проводить по извфетнымь угламъ ф нормали и 
откладывать по нимъ длины ра. а вычислять координаты центровъ 
кривизны соотвьтетвующихь нанесеннымъ на черт № точкамъ мери- 
лана по фориуламь: 

д-р: 608 
В=г— риф. 


Нанеся эти точки и получимъ эволюту, соединяя же каждую 
изъ нихь съ соотеётствующею ей тозкор мернмана получнмЪ нор- 


— М4 — 


мали й длины р,. ЧТо послужить контролемъ какъ вёрности Черкел:а, 
такъ и отеутстыя чувствительныхь оптибокь въ вычисленит. 

Примфромъ такого чертежа служить приложенный къ этой статьЪ, 
который быль исполненъ, принявъ а = 200 мм., для статьи же умень- 
шенъ фотографически, такъ что @— 177, 8 мм. 


— | 


Фяг, 2. 
Величины @ для разныхь жидкостей стёдующйя: 


Ртуть 
Вода... - 
Керосивъ. . 


Спирть 


Формы капилаярныхь поверхностей приведены въ снимкахъ съ 
тицательно исполненныхь чертежей, между прочимъ въ т Г „Рорщаг 
1есёигез ап адтеззез“ МУ. Тпотазотга (Тот@ Кеут). Нашь чертежь 
соотвфтетвуеть чертежу 27-ому Томсона. 

Замфтимь также. что формы капель вычислены съ большою 10т- 
ностью и для пихь составлены обширныя таблицы въ сочинении 
Вазио, Ан айетрЕ Ю 5 фе Феочев 0 сарШату псйон. Сатбце. 1985, 
причемь, метода интегрпроватя уравненя’ разработана, знаменитымь 
астрономомь Аламеомъ. 

Метода Адамса аналогична методЪ Шиёрмера и у него на стр. 18 
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приведена формула, отличающаяся лишь обозначешями оть нашей 
основной формулы (13) а именно: 


251 


1 5 Е 
—%=щ («+ А, А 32% +26 ть А -- 


ТДЪ ш есть промежутокъ между посльдовательныхи значешями пере- 
виной независимой # уравнеше предложенное для интегрировавя 
онъ пишеть 


4, 
== 


и черезъ 4.,А4....А 4, обозначены частныя значетя величины 4 и 
ея разностей до 9-аго порядка расположенныя въ нижней косой строкь 
таблицы подобной таб. П $ 9. 

Но Адамеъ для вычиелешя вмЪсто этой формулы береть другую 
представляющую разность у/— у; въ этой формулЪ, козффищенты 
при 424,.... АЗуузначительно меныше нежели въ предыдущей волЪлд- 
етв]е чего вычислеше можеть быть произведено съ большею точностью. 
но самый процесеъ сложнЪе. 

Уравнен!е опредъляющее форму поверхности Адамсъ пишеть въ 
такомъ видЪ: 


т.е. выражаеть величины р, и = или при нашемь обозначеши г 
въ доляхъ ращуса кривизны въ вершинъ т.-е. 5, и за начало коорди- 
вать приннмаеть вершину. 

Такимь образоэгь въ табдицахь Вазн(ог'’а нашъ случай соотв®т- 
ствуеть значеню 8-=2 н наше 2 получится придавая 1 кь показан- 
ному въ таблиц. Сличене покажеть что разница не превышаеть 
2 или 3 единицы четвертаго знака (таблицы Ваз огВ“а вызиелены на 
3 знаковъ) такъ что полученный нами результать даже лучше ч®мь 
можно быле ожидать при вычислени на 4 знака и притомъ линейкой. 

$ 2. Скажемъ еще нЪеколько еловъ по поводу тЬхъ многихь урав- 
НОВИЙ математической физики, интегрироване которыхъ приводитея 
въ нахождению такъ назывземыхь фундаментальныее функшИ. Эти 
функщи опредфляются въ каждомъ вопросё обниновеннымъ линей- 
нымъ дифференнальнымь уравнешемъ, заклюзающимъ нейзЗвБФСТНЫЙ 
параметрь, значешя котораго даются нфкоторымъ трансцендентнымь 
уравненемь, подучаемымъ на основани граничныхь условйй. 
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Такъ напримЪфръ, изсаБдоване колебашй упругаго стержня цере- 
минной площади поперечнаго сВченя приводится къ розысканю 
функшИ Х опредъляемыхь уравнешемъ: 


В 
ох] чФх-о (.... “ 


и граннчными уеловями. 

ЗдБеь Хесть н®который параметръ, значешя котораго также тре- 
буется найти, р(т) и 9(2} заданныя функщи. 

Траничныя условя выражаются различно смотря по тому, закрп- 
ленъ, подперть или свободенъ разсматриваемый конець стержня, для 
котораго абсциеса х=0 или х=Ь 

Пусть этоть конешь есть тоть, для котораго х=0, тогда гранич- 
ныя условя выразятся такъ: 


. р ах 
а) если конець закрюялена: при 2-0 должно быть Х--би-т-:-0 
5) „ „  побиерт: „я ”„ 


8) „ „ свободенв: „ „ „ 


Совершенно подобно выражаются эти условйя и для другото конца, 
для котораго 5-1. 

Положимь напримЪрь, что оба конца стержня свободны (корабль 
плавзюц!! на водВ). тогда граничныя условя будуть: 

При х—0 должно быть: 


Хх #Х 

аз = ок и =) 
При 2-1! должно быть: 

2х фХх 


я -=9 и Я =° . а (0% 


Такь какъ уравнеше служащее для опредфлевя Х линейное, то 
общйй его интегралъ есть линейнал функшя произвольныхь постоян- 
ныть т.е. 


Х-=АУВИНСС-ЬГ 
гдъ У, А. СЕ. К частныя рышеня уравненя (*) А, В. С.Б произволь- 


ныя постоянныя. Очевилно, что Г, 2, Е, Г суть функцш какь пере- 
иЪиной х такъ и параметра ^, 
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Граничныя условы будуть, обозначая значками производных по 2: 


АУ" (ол ВИ (о -- Сол -1- Бо = в 
дУ"(ох)-- ВР" (ол) Обо) БР" (Ху ==0 
АРА ВР" Об" БУРЫ) — 0 
АТ"ил-- в2" ах) Об") рт" 0 


Чтобы эти уравневя для 4, В, С, Р допуекали рыненя отлич- 
ныя оть нуля, необходимо чтобы опредфлитель: 


д0) = | Ух, "ох, Рам, У" =о. .. 9. 


Это и есть то транецендентное уравнене, о которомь сказано выше, 
‚и которымъ опредфляются значешя параметра ^. 

Но ясно. что лишь для весьма немногихь частныхъ видовъ функ- 
ЩИ р (2) и 4(2) вопросъ можеть быть проведенъ аналитнчески до конца, 
въ остальвыхь же случаяхь, особенно когда фувкщи р(2) и 4(®} 
заданы графически пли таблично, не только рыпить, но и составить 
въ развернутомъ видЪ уравнешя ([) невозможно и ‚слЬдозательно, надо 
прибфгать къ приближеннымь способамъ, 

Для того, чтобы найти значевя » удовлетворяюцие уравненю {1} 
ять надобности имфть аналитическое выражен!е функши А (А), доста- 
точно имЪфть возможность вычислить частное ея значеше соотвфтетвую- 
щее любому избранному частному значению параметра »№---#. 

Въ самомь дётЬ, если на основаши какихь-либо соображенй, 
изложене которыхъ здЪсь опускаемт, было бы извЪстно, что прибли- 
женное значене \ есть », то вычиеливь А(*), полатаемъ затЬмь 
А=и-- А и вычисляемь А(»--®); если эти величины окажутся раз- 
ныхь знаковъ, то искомое значене № лежить между м и я--®, посль 
чего обычнымь премомъ и найдется сь любой степенью точности. 
Нетрудно видфть, какъ поступать, если бы А (”) и А(я--®) оказались 
одного знака. 

Ходь этого вычиеленя численнаго значешя ДА {7) и вало 
показать. 

Ностоянныя произвольныя можно воображать выбранными такъ, 
чтобы онф представляли частныя значешя фуныши Х и первыть 
трехь ея, производныхь при х-- 0, 

Дая этого стоить только вообразить, что У, Ё, &. Г суть тавя 
частныя рёшешя уравненя 


9 [269 |--ваб9. 


— 118 — 


что при х=0 значешя этихъ функц и ихъ производвыхь таковы: 


ие "о | 


7-0 Р-а 2-0 250 


0 
"о 


Й=.0 И -—0 та. | 


Тань какь Х должно удовлетворять уеловямъ (**) то при тавомъ 
выборЪ функши У, 2, 0. У будеть 


@=си2=0 


и останется только найти функщи Ги #. 
\иобы найти функию Г положим 


яв =1 (2) 
в замфнимь уравнене, коимъ опредЪляетел У снетемою: 
СУ я | 
Ч 14. ] ® 
Начальныя значеня будуть: при х==0 должно быть: 
ен. @ 


Кь системь (3) при усломяхь (4) непосредственно примфнима 
метода Штёрмера, по которой и составимъ таблицу значений! функшй 
Уну на протлженш оть х=0 до х--1, при чемъ попутно получится 
и таблица значений У”. 

Совершеняо также составимь и таблицу значенй 'рункши Я и 
попутно въ нею Я”. 

Некомая фунишя Х будеть: 


ХЕАРТВЕ 
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пользуясь составленными таблицами значешй функщй У” и # нахо- 
димь величины ихъ соотвзтетвуюния 1 =, пусть эти велачины будуть 
Ри 4. Пользуясь значешями У" и " п извъетными формулами вы- 
ражаницимн зависимость между разностями н производкыми состав- 
ляемъ значення У” и Й" при х ==! (если функши р (2) и 942} заланы 
аналитически л це табтично. то эти знатеня найдемь и изъ пред- 
ложеннаго уравнешя продиииеренцирававь его), пусть они будуть 
Ми \, тогда па основанит условий (*°) и (***) должно быть, 


АРВ 0 
АМУВЫ: 


и чваетвое эначене опредълителя А (*) есть 


— 0. 


А(и)= РА 


Вели оно окажется равно нулю, 10 Х==и есть одно изъ нетин- 
пыхь значешй параметра Х, если же нётъ, то надо повторить процессъ, 
положивь ^=я--Ёит дл, примфняя извфстную методу разысканя 
значеня корня близкаго нь данному чиелу и. 

ТакяуЪ образомъ, изложенная метода лаеть практическую воз- 
мольтость при рёшенйг физичеекихь вопросовъ на самомь ЛЪЛЬ со- 
ставлять таблицы и вычиелять значеня фундаментальныхь функц, 
существоваше и свойства которыхъ устанавливаются другими ©посо- 
бами. Примфра численнаго вычислешя мы ве приводимъ, такь какъ 
оно ничЬмъ не огличается оть показанныхь выше и не можеть пред- 
ставить никакихь затруднени, кромф нъкоторой утомительноети по 
своему однообразно. 


Москва. ТАПОГРАЗЧЯ хг ая икЧАТИЯ”, Больш. ГаВаЫ, 1% 1918. 


ТАБЛИЦА УШ метода (Рунге). 
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О ИИ = — — Нк | ито и. 
О Г, И, 134) а ЕЕ и 3 ато ие. и н,_1 ев ® и 5 у у \у 
ны | №" о № Де 5 ” т ОНИ мы НИЕ ее 
| | И ОВ 1 ` 51055 к #11538 Го 03723 в | 9039 ИВ ПавО 
0.100 ' 0.03090 _ 9,49201 _ 0,0615 < ] | 058 | | 008750 Й 19365 ` 0.39906 т Г 56355 | 051058 [0-98 | ооо змова И Оо С в м . Г 2 +10 | 0ю 5685 
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